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Abstract 

| The Gardner method, traditionally used to generate conservation laws of in- 

tegral) le equations, is generalized to generate symmetries. The method is demon- 
strated for the KdV, Camassa-Holm and Sine-Gordon equations. The method 
involves identifying a symmetry which depends upon a parameter; expansion of 
this symmetry in a (formal) power series in the parameter then gives the usual 
infinite hierarchy of symmetries. We show that the obtained symmetries com- 
0^ | mute, discuss the relation of the Gardner method with Lenard recursion (both for 

generating symmetries and conservation laws), and also the connection between 
■ the symmetries of continuous integrable equations and their discrete analogs. 
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1 Introduction 

• — ■ 

More than forty years ago Miura discovered the so-called Miura map [18] . In a subse- 
quent paper Miura, Gardner and Kruskal showed how to use the Miura map for the 
construction of an infinite hierarchy of conservation laws for the Korteweg-de Vries 
(KdV) equation [TH]. For brevity we call this technique the Gardner method. Other 
techniques for the construction of conservation laws have appeared since then, for ex- 
ample Lenard recursion [91 [151 EH El E2] , the Gelfand-Dickey method [3 [101 EZ| and 
the symmetry method [23]. Nevertheless the Gardner method is considered to be the 
simplest and most elegant method for the construction of conservation laws. It has 
been applied to other integrable equations, for instance to the Camassa-Holm(CH) 
equation [7]. 

An analog of the Gardner method for discrete equations appeared recently [251 EE] . 
An infinite number of conservation laws were constructed for discrete KdV (dKdV) 
and all the ABS equations [2] using Backlund transformations (BTs — note the Miura 
map is one of the denning equations for the BT of KdV). Moreover in [25] it was shown 
that it is also possible to use the Gardner method to construct the infinite hierarchy 
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of symmetries for dKdV. This result raises the question whether it is possible to do 
the same for continuous equations. 

The theory of symmetries for continuous equations is well developed. There is 
a direct method which allows computation of all symmetries of a given order for a 
given equation. There are numerous methods to generate the infinite hierarchies of 
symmetries of integrable equations such as KdV, for example, the mastersymmetry 
method [8], Lenard recursion j22j E2 E3] and Lax operator methods [10j HSl [31]. The 
method we will present in this paper has much in common with the method of the 
resolvent, described, for KdV, in section 3.7 of [5]. 

The goal of this paper is to present the Gardner method for symmetries of contin- 
uous partial differential equations (PDE). It is a little more subtle than the Gardner 
method for conservation laws, but it is still simple and elegant. The structure of this 
paper is as follows: In section 2 the Gardner method for symmetries is presented for 
KdV, CH and sine-Gordon (SG) equations. In section 3 we prove that the obtained 
symmetries commute. In section 4 we show the connection of the Gardner method and 
Lenard recursion. In section 5 the connection between the symmetries of continuous 
and discrete equations is described. Section 6 contains some concluding comments and 
questions for further study. 



2 Gardner method for symmetries 
2.1 KdV 

The potential KdV (pKdV) and KdV equations have the form 

u t - \u\ - \u xxx = , (1) 

4>t - 34>4>x - \4>xxx = o (2) 

related by 4> = u x . The standard Backlund transformation for u is u — > u a = u + v a 
where v a satisfies 

v a ,x = a - 2<p - v 2 a , (3) 
v a ,t = ~\4>xx + {a + 4>){a - 2(j) - v 2 a ) + <f> x v a . (4) 

The index a indicates that v a is the solution of (J3]) with parameter a. The system ([3])- 
(j3]) is consistent if <fi satisfies the KdV equation ([2]) and has a one-parameter family of 
solutions. By v ai , v a2 we denote different solutions of ©-dH) with the same parameter 
a. There is an algebraic way to describe the action of repeated Backlund transforma- 
tions with different parameters [6]. This relation is called the nonlinear superposition 
principle and has the form 

U a ,P = U^ . (5) 

Up - U a 

Here u a = u + v a is obtained by application of the BT with parameter a, up = u + v@ 
is from application of the BT with parameter /3, and u at p is the result of application 
of both, in either order, as they commute. Writing (3 = e + a, relation §5§ can be 
rewritten 

e e 
w QiQ+e = iiH = «H — . (6) 
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For small e we can expand v a+e as a power series in e. To lowest order, v a+€ satisfies 
the same system of differential equations as v a . But this does not mean that to lowest 
order it is the same as v a , since, as we have explained, the BT has a one-parameter 
family of solutions. Relation (E]) can thus be rewritten in the form 

u a , a+e = u + ^ + 0(e 2 ) . (7) 

where v ai and 

V0L2 £1X6 distinct solutions of ([3j)- ([4j) for the same parameter value a. 
Relation Q describes an infinitesimal continuous transformation of u which is also 
a solution of pKdV. In other words, this is a Lie symmetry for pKdV. We denote it 
X(a) = Q{<y)-§^ where 

Q(a) = i . (8) 

We call X(a) the general symmetry. Note that while the explanation we have given 
above for why X(a) is a Lie symmetry is perfectly rigorous, there is a more direct 
technical proof: If Q is defined as in §E§ and v ai and v a2 are both solutions of (J2])-® 
then it is a technical exercise to check that Q satisfies the equation 

Qt - - \Q XXX = , (9) 

this being the defining equation for infinitesimal symmetries of pKdV ([!]) . 

The next thing to do is to observe that if we could solve ©-((1]) to write as a 
function of <fi (or u) and (a finite number of) its derivatives, as well as the constant of 
integration, then we could rewrite (jBJ) in terms of <fi. This cannot be done explicitly, 
but it is possible to write a formal solution of (]3])-® to give v a in terms of as a 
formal series in decreasing powers of a 1 / 2 . Explicitly we have 



1/2 V Vx Vxx + ±y <±>xxx T Vxxxx T "V ' T *-*H"+>xx ^ / _3\ 

a ~aV* + 2a 4a 3 / 2 8rf +U V* > 

(10) 



2 



or 



oo 

Ci 



i Oil 

4 = 1 " 



where 



n-l 



= TT , C n+1 = ( Cn,x + ^ CjCn-j \ , n = 2,3... . (12) 



2 ' 2 



8=1 



At first glance this seems insufficient for our purpose, as it only gives a single solution 
of equations dm)-®, with the prescribed asymptotic behavior for large \a\. But there 
evidently is a second solution in which a 1 / 2 is replaced by —a 1 / 2 . Thus we can take 



v ai = a 2 



^+£4, *, = -a*+f;fc^, (i3) 



1=1 " i=l 



with the c, defined as before to obtain 



<X«) = 2^1/2/ -1 ~ x ( 14 ) 



(i + e; 



OO C2j + 1 
+ 1 



1 (fixx 60 (pxxxx 

+ 4O0 3 + 1O0 2 , + 20<jxpxx 



2a 1 / 2 V a 4a 2 16a 3 
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This expansion gives the infinite hierarchy of symmetries of pKdV. The first few sym- 
metries take the form 



X - 8 

OU 

Y 9 
Al = U x — , 

OU 

d d 
X 2 = (u xxx + 6u x ) — = Au t — , 

d 

X 3 = (u xxxxx + 40u x + I0u 2 xx + 20u x u xxx ) — . 
The corresponding symmetry for KdV is 

Y(a) = • (15) 

(Note that in general it is nontrivial that a local symmetry of pKdV should give a 
local symmetry of KdV, as u is nonlocal in 0, but since all the symmetries for pKdV 
we are considering are determined by <fi this is not an issue here.) Using (jBl),© and 
(fT3"|) we have 



V^OO C2j+2 

(v ai -v a2 y v ai - v a2 a 6 ' 2 ( 1 _|_ ^ \ 

Expanding in powers of -, or just differentiating the relevant formulas for pKdV, gives 
the first few symmetries for KdV 



Y-i 



Y 2 = + 12<p<f) x ) — = 4<f) t — 



d 

Y 3 = (4>xxxxx + 12O0 2 X + A0(f) x (j) xx + 20(j)(J) xxx ) — . 

Thus we can generate explicit formulas for the infinite hierarchy of symmetries of 
pKdV or KdV using ffl2|) to find the c$ and then expanding (I14p for pKdV or (TT6]) for 
KdV in inverse powers of a. 

In fact the function Q(a) that has emerged as a generating function for symmetries 
of KdV can be identified with the "resolvent" of KdV (see [5] section 3.7). Writing 
v a = equation (jSJ) becomes the Schrodinger equation 

Ipaxx = (a- 20) ^ Q • 

Using dHJ) and the fact that the Wronskian of two solutions of the Schrodinger equation 
is a constant, we find that Q(a) can be identified with the product of two solutions 
of the Schrodinger equation, which is the resolvent. However, the proof we have given 
above that Q(a) is a symmetry of pKdV is new, and as we shall see, this new proof 
allows generalization to other equations, as well as a simpler proof of properties such 
as commutativity. 
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2.2 CH 

In [7] the conserved quantities of CH were derived from those of the associated Camassa- 
Holm (ACH), introduced in [23] and studied further in [T21 E2 El M\- The ACH 
equation has the form 

P ( Pt\ P 2 



* = f = i(j). t (17) 

In [29] a BT for ACH was given in the form p — >■ p — s x where s satisfies 

s 2 a , . 

s x = + - + P, 18 

pa p 

s t = - s 2 + —s + a(a-2f). (19) 
P 

The superposition principle for ACH is 

aj^M^sl) . (20) 

V ps a — asp J x 

Here p a = p — s a x is obtained by application of the BT with parameter a, pp = p — sp. x 
is from application of the BT with parameter /3, and p a> p is the result of application 
of both. Writing j3 = e + a, and doing the expansion of (1201) around e = we obtain 



( ° 2 SaiSa2 ) 



s ai s 



"2 / X 



where s ai and s a2 are distinct solutions of (I18|I19I) for the same parameter value a. 
Relation (j21[) describes an infinitesimal continuous transformation of p which is also a 
solution of ACH. We denote it X(a) = Q(a)-^ where 

Q(a) = f a2 ~ s ^ Sa A = K^L±^al. (22) 

Using ( TT7|ll8lll9p it can be checked that the coefficient of Jj in the prolongation of 
X(a) is 

Q = a ( Sai + Sa2 - Pfx ^ . (23) 

Furthermore Q,Q satisfy the equation for infinitesimal symmetries of ACH 

Qt = 2pf x Q+p 2 Q x . (24) 

To generate the symmetries we apply the same idea as in the previous section. 
Namely, we write a formal solutions of ( Fl8|T9|) as a series in increasing powers a 2 



00 

En 

n=l 



where 



P 

2 ' ~'^ x 2 ' 2p 



n-2 



i=0 
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The second solution of ( 1181191) can be obtained by replacing a 2 by — a 2. Thus we take 

00 00 
s ai = ^s n a% , s a2 = ^ s n (~a^) n . 



n=l n=l 

Plugging this into (j22p we obtain 



(25) 



The expansion of (125]) around a = gives the infinite hierarchy of symmetries of ACH. 
The first few take the form 

Xi = PxS- , 
op 

X 2 = ( 2p xxx - 12§ - 6^ + 3^) I- , 

\ £r p p A J Op 

~Y ( ^ cyPxPxxxx ^PxxPxxx ^-^-PxPxx ^PxPxxx ^Pxxx 

^3 I "zPxxxxx ^ 1 



5 p p 

PxPxx (28-19pl) | 3^(16-40^ + 9^ 9 
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p3 4^4 y ^? 

As far as we are aware these symmetries are new. 

2.3 SG 

The equation SG 

u xy = sin(u) 

can be brought to rational form by the change of variable u = 2i ln(z), giving 

zzxy - z x z y = -(z 4 - 1). (26) 

The BT for (127)]) jl] is 2 — > z Q where z a satisfies 

(z a z)x = ^(zl- z 2 ) , zz a , y - z y z a = ^-(z 2 z 2 a -l). (27) 

The corresponding superposition principle is [UJ 

^ = s-(a + fl gfe— 33 (28) 
pz a — az/3 

Plugging j3 = —a + e and expanding (1201 around e = we obtain 

Z{Z ai + <Z a2 ) 2 \ /on\ 

z Q = z + e— - + 0(e), (29) 

where z Q1 = z a and z Q2 = — z_ a are distinct solutions of (1271) for the same parameter 
value a. Relation ([2"9l describes an infinitesimal continuous transformation of z which 



is also a solution of SG. We denote it X(a) = Q(a)S- where 



= z(z ai +z a2 ) 



One can check that Q satisfies the equation for infinitesimal symmetries of SG 

Qz X y ~t~ xy QxZy ^xQy Z Q 0. 

We write a formal solution of (1271) as a series in decreasing powers of a 

oo 

EVn_ 

where 



a' 

n=0 



_i_ 1 ST^n 

. . _ Z x V n ~T ZV nx Z^i=l v i v n+l— i _ _ 

v = ±z, v 1 = ±2zz x , v n+1 = , n = 2,3.... 

There are two distinct expansions of z a . They can be used as two different solutions 
z ai and z a2 of (1271) . Plugging z ai and z a2 into (130]) and expanding around a = we 
obtain an infinite hierarchy of symmetries. The first few symmetries take the form 

X z 8 

-<M — z x — ; 

OZ 

V f r,ZxZ X x\ d 

~ \ Zxxx ~ — ~ — ) ~Q~ ' 

2,0z x z xx \0z x z xxx ^-0z xx z xxx \0z x z xx &z x z xxxx \ d 

3 I Z xxxxx + - "T 



z 2 3 z J dz 

These symmetries coincide with ones obtained in [22]. However, as far as we know 
the derivation above is new and represents a substantial simplification over existing 
techniques to generate the hierarchy. It is straightforward to verify that the j'th 
symmetry for all equations (in the above numbering) depends on x-derivatives of u, p, z 
up to order 2j — 1, and that the coefficient of the highest derivative is nonzero, thus 
guaranteeing nontriviality. 



3 Algebra of the symmetries 

In this section we determine the algebra of the symmetries obtained in the previous 
section. We prove that symmetries generated by the Gardner method for KdV, CH 
and SG commute. In the previous section we showed that symmetries for KdV can be 
obtained by the expansion of the general symmetry X(a) in a series around a = oo. 
In order to prove that these symmetries commute it is enough to show the general 
symmetry commutes with itself but with a different parameter, namely that 

[X{a),X{P)] = 0. (31) 

(We refer in this section only to the symmetry for pKdV, the result follows for the 
general symmetry of KdV Y(a) by a simple prolongation argument.) Here 

Id Id 

*(<*) = 7. —^7. X{J3) = Z -75T- ( 32 ) 



v, 



a 2 



du ' vp 1 — V/3 2 du 



To compute the commutator we need to know how X(a) acts on v@ 1 and Vp 2 and how 
X((3) acts on v ai and v a2 . This is currently not clear. However we can also write 
X(a),X(P) in the form 

Id Id 

X(a) = -f , X(J3) = -f . (33) 

u ai — u a2 du u$ x — up 2 du 
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The relevant prolongations of X(a) and X(f3) are then 

X{a) = 4- + tt— + 7T— , 



1 


d 


1 

f- 


d 


U ai fa — 


U a2 fa ® U fa 






1 


d 


1 


d 


u aifa — 


^oi/32 du ai 


h 

U 0.2 fa 





Ufa — Ufa ou u ai/3l — u ai p 2 ou ai u a2/3l — 
To determine the commutator of X(a) and X(f3) we compute 



X{a)Q{0) - X{p)Q{a) = X(a) I - X{$) 

\u Pl - Ufa J 

1 / 1 



{Ufa — Ufa) 2 \ U ai fa — U a2/3l U ail32 — U a2/ 3 2 

1 / 1 1 

+ 



(36) 



(u ai U a2 ) y U ai fa U ai fa U^fa U a2 fa 

In order to simplify this expression we have to use relation (J5J). Since we have two 
different BTs of u for parameter a and two different BTs of u for parameter (5 we need 
relation (|5]) in the four different forms 

a-P , 

u a.ifa — u-\ , i,j E {1,2} . (37) 

u ai - Ufa 

With the help of ( 157)) the right hand side of ( 155)) simplifies to zero, thereby establishing 
(l3Tj) . Thus the symmetries obtained by the Gardner method for pKdV (and KdV) 
generate an abelian group. 

Let us prove a similar result for ACH. For this we need to compute the prolongation 
of the infinitesimal generator of the general symmetry. This can be done by the change 
of variable 

u = pdx. (38) 



The superposition principle for ACH can be rewritten in the form 

f(a-/3)(af3- {u-u a ){u -up)) \ , . 

U a ,0 — U — — ; r r . {6U) 

\ p{u — u a ) — a{u — up) J 

This expression is equivalent to quad-graph equation Qls=i in the ABS classification 
[2]. Using (138)) . the general symmetry of ACH is 

a 2 — (u — u ai )(u — u a2 ) d 



X (a) 



du 



In order to prove that the symmetries commute it is enough to show that ( 15T|) holds 
for the general symmetry of ACH. The relevant prolongations of X(a) and X{/3) are 

- a 2 - [u - u ai )(u - u a2 ) d a 2 - (Ufa - Ufa ai )(ufa - Ufa a2 ) d 
X{a) = — h 



+ 

X(J3) = 
+ 



du Ufa a2 - up iai dufa 

® 2 - ( u fa ~ up 2ai )(ufa - Ufa a2 ) d 

U faoi2 U faa 1 dup 2 ' 

/3 2 - (u - Ufa)(u - up 2 ) d f3 2 - (u ai - Ufa ai )(u ai - Ufa ai ) d 



up 2 - Ufa du Ufa ai - Ufa ai du ai 

f3 2 - (U a2 - Ufa a2 )(u a2 - Ufa a2 ) d 



Ufaa 2 — Ufaa 2 du 



a 2 
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Relation f l39|) can be presented in four different forms which connect different BTs for 
CH 

/ (oj - -{u- u ai )(u - Uf } )) \ . 

V /3{u- u ai ) - a{u- u^) J 

With the help of this we obtain that (1311) is true for the general symmetry of ACH. 
Thus the symmetries obtained by the Gardner method for ACH generate an abelian 
group. 

The proof that symmetries obtained by the Gardner method for SG commute is 
very similar to the above. The prolongations of the general symmetries X(a) and 
X((3) for SG are 

\ z(z ai + Z a2 ) 8 Zpi\Zpia\ + Zj9ia 2 ) 9 Zp 2 {zp 2ai + Zp 2(3t2 ) 8 

■X-(a) = —, r-— H p r— — h 



^ai Z a2 

lai z 0ia-2 ) dz f32 



Relations which connect different BTs for SG can be obtained from (1281) 



^ A = g -(a + fl fe ^ i,je{l,2}. (41) 

This helps us to obtain that f l3~Tj) is true for the general symmetry of SG. Thus the 
symmetries obtained by the Gardner method for SG also generate an abelian group. 



4 Connection of Gardner method with Lenard re- 
cursion 

As we mentioned in the introduction there are several methods for generation of sym- 
metries for PDE. One of them is Lenard recursion [221 1231 [21] • We intend to show 
that Lenard recursion gives the result equivalent to the Gardner method. We present 
a brief explanation of Lenard recursion for the KdV equation. 

The Lenard recursion is based on the fact that KdV (j2J) is a bi-Hamiltonian system. 
Namely KdV can be presented in two forms: 

u t = (42) 



where H and Hi are 



/>3 



and 



#o = / ( y - f ) dx, E x = I ^dx, 



P - — P - -— + 20— + 

dx ' 4 dx 3 dx 



The expression ^ denotes the variational derivative. If g — J Gdx is a conserved 

quantity for KdV equation then both Q = -Pof|) Qi — Pi| are the characteristics of 
symmetries. 



9 



Evidentiary 

PiPo Qo — Qi- 

Thus the 'recursion operator' for symmetries 



R = P1P0 1 

maps Hamiltonian symmetries to Hamiltonian symmetries, and can be used to generate 
the infinite hierarchy of symmetries of KdV. Though it is necessary to check that at 
each step the operator P _1 can be applied. 

In section 2 we showed that KdV has a symmetry with characteristic 

1 \ v ai + v a2 



Q(ct)x = 

Applying R to Q(a) x we obtain 

RQ(a) x = P^ 1 ( 1 ) = Pi \ = a ^+ v ^ = a Q{a) x . (44) 
Writing 



a 1 / 2 *■ — ' a ' 

i=0 



and substituting into (144}) we obtain 



R I Qi \ _ Qi 

\i=0 / i=0 

Comparing coefficients of powers of a we get 

RQi = Q i+1 , 2 = 0,1,2,.... 

Thus Lenard recursion gives the same symmetries as the Gardner method for KdV. 
Further more, in our approach we see immediately that all the characteristics are x 
derivatives. 

We briefly demonstrate that Lenard recursion for conservation laws can be derived 
from the Gardner method. In the Gardner method the infinite hierarchy of conserva- 
tion laws for KdV is obtained by expanding the conservation law 

d t (v ai - v a2 ) + d x (-(u + 6)(v ai - v ai )) = 0. (45) 

into a series in a with the help of f[T51) [T§] . The conserved quantity associated with 
(145]) has the form g = J (v ai — v a2 )dx. We now prove that the variational derivative 
of g is 

5q 4 , . 

The proof for ( 146]) is as follows. We know from ([3]) that 



v auX + v a ^ x -2a + v 2 ai + v 2 



a 2 
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The variation of this expression is 

Sv ai)X + 5v a2)X + 2v ai Sv ai + 2v a2 5v a2 = -A8cj). 

Dividing this by v ai — v a2 we obtain 

2v 

| L til CKl ' UC'2 UL'l ^ 



8v aux + 8v a3>x 2v ai Sv ai + 2v a2 5v a2 



This expression can be rewritten differently as 

(Sv ai + 5v a2 )(v ai + v a2 ) ( Sv ai + Sv a2 \ 2v ai 5v ai + 2v a2 8v a2 



_|_ Q I " v «l 1 " ""2 j _|_ ""«r"«l 1 -"«2" v "2 i t 

After simplification we obtain 

0«ai - = 4 \-Ox[- 

^«2 V 

Thus 

4 



5 / (u ai — v a2 )dx = / 

J J "(Jl ^«2 

and ( |46|) is proved. We have already seen in (144]) that 

Pi = 0iQ{a) x = aP - 



Thus 

d ^ p 5 9 

or 

o S 9 $9 

where R = Pq X P\ is the recursion operator for (variational derivatives of) conserved 
quantities of KdV. 

5 Connection with symmetries of difference equa- 
tions. 

Relation (J5J) 

a — 

u af3 = u + 



u a - Up 

connects solutions of the pKdV equation. It can be embedded on the quad-graph, 
the planar graph with quadrilateral faces [2j [3j [201 EQ]. The embedded equation is 
called the discrete Korteweg-de Vries equation (dKdV) and has the following form 

01 ~ P l A n s 

"1,1 = "o,o H • (47) 

"i,o - "0,1 

Here k, I G Z 2 are independent variables and u = u(k, I) is a dependent variable that 
is defined on the domain Z 2 . We denote the values of this variable on other points by 
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Uij = u(k + i,l + j) = S l k Sjuo t o, where Sfc, Si are unit forward shift operators in k 
and I respectively. Let us explain the connection between fH7|) and ([5]). Shifts U\$ and 
uo,i are BTs of u = ito,o for two different parameters a, (3. In general Uij denotes the 
application of the BT i + j times to u = u 0j o, i times with parameter a and j times 
with parameter (5. 

Properties of dKdV are closely related to properties of pKdV. In particular, since 
the continuous general symmetry of pKdV discussed in the previous sections is an in- 
finitesimal version of a Backlund transformation, and Backlund transformations com- 
mute, the general symmetry must also give a continuous symmetry of dKdV, and we 
wish to identify its generator. 

So far we have identified the functions u ai , u a2 in the pKdV general symmetry 

x ^ - < 48 > 

as different Backlund transformations of u with the same parameter a. In fact we can 
just as well identify one of them as a forward Backlund transformation and another 
as a reverse Backlund transformation. To demonstrate this, observe that a reverse 
Backlund transformation of u is a function u such that u = u + v where v solves the 
system 

v x = a - 2(j) - v 2 , (49) 
v t = -\4>xx + (a + $)(a-24>- v 2 ) + 4> x v , (50) 

and (p = u x . Replacing in fT4T?l) - floll by — v x , and using the first equation in the 
second to eliminate all x derivatives of v, the above system becomes 

— v x = a — 2cj) — v 2 , (51) 
-v t = -\(f)xx + {a + 4>){a - 2(f) - v 2 ) - <j) x v . (52) 

Comparing with ©-(jl]), we can identify v with — v a , where v a is the function appearing 
in the forward Backlund transformation. In other words, forward and reverse Backlund 
transformations are identical. 

In the language of the quad-graph, inverse shifts W-^o and «o,-i are reverse BTs. 
So we can identify u a2 with -u_i i0 and u ai with -ui j0 and the quad-graph version of 
X(a) is thus 

Xt = J- . (53) 

Similarly the quad-graph version of X(j3) is 

X 2 = /- . (54) 

These symmetries were already found in [27]. Furthermore, in [25] it was shown that 
equation fT4Tl) embedded in three dimensions has another symmetry of the form 

X = ^— . (55) 

^0,0,1 — ^0,0,-1 au 0,0,0 

This evidently also has its origins in the general symmetry. In [25] it was shown how, 
by a suitable expansion, this symmetry gives the infinite hierarchy of symmetries for 
the dKdV equation. 
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As we noted before the CH equation after the change of variable u = J pdx is 
related to quad-graph equation Qls=i- Namely, the superposition principle for CH 
f )39|) is equivalent to quad-graph equation Ql$=i in the ABS classification [2] 

6» 2 (iti,i - u 0A )(u 0t0 - ui, ) - - u lt0 )(u Oj0 - u 0j i) + (6»i - 9 2 )9 1 9 2 = 0. 

So from the general symmetry of ACH we can obtain the symmetry of Ql$=i. 

By doing the same computation as we did for pKdV we can identify forward 
Backlund transformation and u a2 as a reverse Backhand transformation. Thus 



X 1 
X 2 



O? - (U fl - Uifi) (m ,o - m-i,o) \ d 



« 2 - (u ,o - Uqa)(Uo o - Uq i) \ d 



u -x - u ,i J du 0fi 



are the symmetries of Qls=i- 

For SG we can identify forward Backlund transformation and 

reverse Backlund transformation. From general symmetry of SG we obtain that 

v ^0,0(^1,0 — z -i,o) 9 

A 1 ■ 



(^1,0 + ^-1,0) ^0,0 

Y- _ Zq,o(Zq,1 — ^0,-1) <9 
(^0,1 + ^0,-1) ^0,0 

are the symmetries for equation 

c^V^o.o^o,! — ^1,0^1,1) — ^(^0,0^1,0 — ^0,1^1,1) — 0. 
This equation is equivalent to H3s=o i n the ABS classification. 



6 Conclusion 

In this paper we have introduced the Gardner method for generation of the infinite 
hierarchy of symmetries of integrable equations, using KdV, CH and SG as examples. 
The method involves identifying the general symmetry X(a) from the superposition 
principle for BTs of the equations studied, followed by a suitable expansion in powers of 
a. The method is both mathematically elegant and computationally efficient. We have 
shown how to use our formalism to prove the symmetries commute, explained the origin 
of Lenard recursion relations, and explored the link with integrable lattice equations. 
The fact that integrable lattice equations arise as the superposition principle for BTs of 
continuum equations is well-known — indeed the celebrated Q4 lattice equation in the 
ABS classification was first written down by Adler [1] as the superposition principle 
for the Krichever-Novikov equation. However, we believe the link between the Ql$=i 
lattice equation and CH given in this paper is new, as are the symmetries derived for 
the ACH equation. 

We expect the method to be applicable to other integrable equations, and it will 
be interesting to see more examples developed. It would be particularly interesting to 
see examples in which infinite hierarchies of symmetries and/or conservation laws can 
be constructed by one method, but not by another (though this may be difficult to 
verify) . 
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